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In a normal-metal/insulator/superconductor (NIS) tunnel junction refrigerator, the normal-metal 
electrons are cooled and the dissipated power heats the superconducting electrode. This paper 
presents a review of the mechanisms by which heat leaves the superconductor and introduces over- 
layer quasiparticle traps for more effective heatsinking. A comprehensive thermal model is presented 
that accounts for the described physics, including the behavior of athermal phonons generated 
by both quasiparticle recombination and trapped quasiparticles. We compare the model to mea- 
surements of a large area (> 400 fxm^) NIS refrigerator with overlayer quasiparticle traps, and 
demonstrate that the model is in good agreement experiment. The refrigerator IV curve at a bath 
temperature of 300 mK is consistent with an electron temperature of 82 mK. However, evidence 
from independent thermometer junctions suggests that the refrigerator junction is creating an ather- 
mal electron distribution whose total excitation energy corresponds to a higher temperature than is 
indicated by the refrigerator IV curve. 

PACS numbers: 74.70.Ad 74.50.+r 07.20.Mc 85.30.Mn 85.25.Na 



I. INTRODUCTION 

In a normal-metal/insulator/superconductor (NIS) 
tunnel junction biased near the superconducting gap en- 
ergy A, the single quasiparticle tunneling current trans- 
fers heat from normal metal electrons to the supercon- 
ductor. This transfer enables refrige rato rs that can cool 
electrons from 300 mK to ~100 mK,'il2l and can cool ar- 
bitrary payloads as well. For example, NIS refrigera- 
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Cooled Electrons 

FIG. 1: Cross-sectional sketch of an NIS refrigerator with 
overlayer trap. A quasiparticle (red dot) is depicted diffusing 
through the superconductor to the overlayer trap. Layers are 
labeled O for the overlayer trap, S for the superconductor, 
for the cooled normal metal layer, Ii and U for insulating 
layers. Because U is fabricated independently of U, it can 
have a different transparency. 



'Contribution of a U.S. government agency, not subject to copy- 
right. 



tors have been used to cool a macroscopic germanium 
thermometer^ and a superconducting transition edge x- 
ray detector.E! The performance of these refrigerators is 
limited, in part, by heating of the superconductor due to 
the dissipated power IV, and the heat removed from the 
normal metal. The impact of this heating on NIS refrig- 
erator performance is often characterized by the fraction 
/3 of the power deposited in the superconductor that re- 
turns to the normal metal as an excess load.^ 

Previous efforts to model the heating of the supercon- 
ductor in NIS refrigerators to predict /3 have included 
quasiparticle diffusion, trapping and recombination. For 
example, Ullom and Fisher numerically solved a differen- 
tial equation for the excess quasiparticle density vs po- 
sition in the superconductor.i^ Rajauria et aP used ap- 
proximations to solve similar differential equations an- 
alytically and introduced a finite quasiparticle trapping 
rate. In that work, a parameter equivalent to /3 was cal- 
culated, and the agreement with experiment was within 
a factor of 3 to 10. 

In this paper, we expand upon previous work by pro- 
viding the most comprehensive model of NIS refrigera- 
tion to date. We add a new form of quasiparticle traps, 
referred to as overlayer traps, to both the model and de- 
vices. We model not only the superconductor quasipar- 
ticle temperature, but also the overlayer trap electron 
temperature and the temperature of the phonons in the 
metal layers that make up the NIS refrigerator. We also 
account for the athermal behavior of excitations with en- 
ergy ^ k\jT created by quasiparticle relaxation. Our 
implementation of this model easily handles changes in 
nearly every input parameter, allowing us to examine a 
large area of parameter space and design the next gener- 
ation of NIS refrigerators. We present measurements on 
a large area NIS refrigerator that agree with the model 
predictions over a large temperature range. The refriger- 
ator IV curves are consistent with cooling from 300 mK 
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FIG. 2: Block diagram representation of the thermal model. 
Roughly, power is deposited in the superconductor and leaves 
by trapping to the overlayer traps. Power leaves the electrons 
of the overlayer trap by coupling to phonons in the metal 
layer, and the phonons couple to the substrate (the bath) via 
a boundary resistance. There is additional complexity due 
to quasiparticle recombination in the superconductor and the 
behavior of athermal phonons created by recombination and 
trapped quasiparticles. 



to 82 mK. As discussed in Sec |VI| the refrigerator may 
be creating an athermal electron distribution so the in- 
terpretation of the cooling results is not straightforward. 



II. MODEL OVERVIEW 

We model a single NIS junction, which makes up half of 
an SINIS refrigerator, as shown in Fig[l] The important 
systems in the device are the cooled normal metal elec- 
trons N, the superconductor quasiparticles S, the over- 
layer electrons O, and the combined metal layer phonons 
P. Roughly speaking, the power flow is as follows: a 
power Ps is deposited in the S layer by the NIS refrig- 
erator junction, which increases the quasiparticle density 
locally above the junction. These quasiparticles may dif- 



fuse and recombine, but the majority are trapped into 
the O layer. At this point, they become electronic exci- 
tations with energy A :s> k\,T. Various processes allow 
these electrons to relax, and the majority of the energy 
couples via electron-phonon coupling to the P system. 
Because all three metal layers are made of Al, or AlMn, 
with only thin oxide layers between them, we assume 
that there is one phonon system shared by all three metal 
layers. Finally, the phonons in the metal layers relax by 
coupling to the substrate phonons. Systems in the model 
are connected by both thermal and athermal processes. 
The model is summarized as a block diagram in Fig [2] 



Section [III] discusses NIS tunneling. Section [IV] de- 
scribes quasiparticle trapping, and a detailed discussion 
of the rest of the physics underlying the model is found 
in Appendices [Xf|C| Table |l] shows all the model param- 
eters and typical values. Section [V] compares the model 
predictions to measurements on an NIS refrigerator de- 
vice. 

The model consists of four coupled equations, three of 
which are position dependent. Equation [T] describes the 
overlayer trap electron temperature Tq. The terms on 
the right hand side from left to right are due to trap- 
ping from iS", coupling to P, and recombination phonons. 
Equation [2] describes the excess quasiparticle density in 
the superconductor rtox- The terms from left to right are 
injection by the junction, trapping to O, recombination, 
and trapping to the side-traps. Equation [3] describes the 
temperature of the phonons in the combined metal layers 
Tp. The terms from left to right are coupling to O, cou- 
pling to N layer, and coupling to the substrate. Equation 
|4] is a power balance equation for the electron tempera- 
ture Tn of the N layer. The terms from left to right 
are due to the tunnel junction (Pn is negative during re- 
frigeration). Joule heating, coupling to P, two-particle 
tunneling (Andreev reflections), recombination phonons, 
phonons generated by trapped quasiparticles, and either 
stray power or power from a payload. 
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The variables to, tg, tp, and tN are the thicknesses of 
the O, S, P, and N layers, kq is t he e lectronic ther- 
mal conductivity of the O layer (Eq |A2[ ), the IIx vari- 
ables describe the relaxation branching ratios of trapped 
quasiparticles (Eq 17 1, the A^-y variables describe the 



probability of abs orpt ion of athermal phonons of energy 
X in layer y (Sec C4), T^{nl^ + 2nex»-th ) is the recom- 
bination rate of excess quasiparticles (Eq |B6[ ), nth is the 
thermal quasiparticle density due to the bath tempera- 
ture, I?s-i is the quasiparticle diffusion constant (Eq B3 1, 
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g and 5side-trap are functions describing the location of 
the refrigerator junction and side-traps, Kp is the metal 
layer phonon thermal conductivity (Eq C2), T^trap is the 
quasiparticle trapping rate (Eq 15l), T'n-n and Pn-n are 
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electron phonon coupling terms (Eq |A3[ ), Vs is the power 
deposited in the superconductor by the NIS junction (Eq 
13), and Pamm is power flow across the acoustic mis- 



match between the Al layers and the substrate (Eq CI I 



Script V terms have units power per unit area or volume, 
where non-script P terms have units of power. All equa- 
tions are solved numerically, and Eqs[T]^have boundary 
conditions dT/dx = at both a; = and x = Xend^ 

In the power balance equation (Eq|4| power Pn is de- 
posited in the N layer by the NIS junction (Pn is neg- 
ative during refrigeration). The temperature reduction 
is limited by an electron-phonon coupling power Pp-o, 
Joule heating I^gRpud, two-particle tunneling dissipation 
l2Vh, incident phonon power from quasiparticle recombi- 
nation Pqp_rocomb = A2A-NiN^NS ArR(n^^ -|- 2ncxnth), 

and incident phonon power from trapped quasiparticles 
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^3A/4-N^trap^NS- Here, Ppad is the 



resistance of the current path through the iV layer, /ns is 
the total tunneling current, Vh is the junction bias volt- 
age, and I2 is the two particle tunneling current (Eq[9|.l^ 
An additional power Pexcess niay be present, such as stray 
RF power or power dissipated by a cooled payload. A 
term /3Ps is often included in Eq |4] and used as an em- 
pirical parameter that accounts for heating of the super- 
conductor. For the purposes of this paper we use (5 = 
and account for heating of the superconductor explicitly. 
One goal of the comprehensive thermal model is to pre- 
dict /3 from measurable NIS design and material param- 
eters. When we compare model and experimental results 
in Section [V] we use /3 as a comparison metric. 

Figure [3| shows the temperatures of the 5, O and P 
layers as calculated by the model for parameters in Ta- 
blejlj Figurefflshows the terms in the normal metal power 
balance (EqpJ) calculated by the model for different val- 
ues of the overlayer trap barrier resistance area product 
7^so• 
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FIG. 3: Temperatures of the superconductor quasiparticles, 
overlayer trap electrons and the phonons of the combined 
metal layers vs position shown from a; = to a:: = 50 /im 
calculated with the comprehensive thermal model and pa- 
rameters from Table |l] The pink bar on the x-axis shows the 
location of the refrigerator junction, and therefore, of quasi- 
particle injection into the superconductor. 



iV(0) is the two-spin density of states at the Fermi en- 
ergy in the same material in the normal state, E is the 
energy relative to the Fermi energy, A is the BCS energy 
gap, and 7 is a unitless factor that describe deviations 
from the ideal BCS superconducting density of states. 
The parameter 7 introduces states below the gap energy 
A and is used to account for observed sub-gap tunneling 
currents greater those predicted with the ideal BCS den- 
sity of states (7 = 0). Pekola et aP suggest environment- 
assisted tunneling as the mechanism responsible for finite 

7- 

The quasiparticle density n in a superconductor at 
temperature Tx is 



III. SUPERCONDUCTORS AND TUNNELING 



A. Quasiparticle density and density of states 



n(rx) = 2iV(0) / UE)iy{E%=odE 



(6) 



This section will briefly review a subset of proper- 
ties of superconductors which are vital for understanding 
their behavior in NIS refrigerators. In a superconductor 
at zero temperature, all the conduction electrons form 
Cooper pairs which can carry electrical current with zero 
resistance. At finite temperature, some of the Cooper 
pairs are broken and each broken Cooper pair yields two 
quasiparticles. The effective density of states of quasi- 
particles is given by N{0)h'{E), where 



i.{E) = 



Rc 



/ E/A - 



(5) 



where fx{E) ~ (^^E/kbT^ ^ ^j-^g Fermi function at 

temperature Tx, and fcb is Boltzmann's constant. To ac- 
count for equilibrium behavior, we write the total quasi- 
particle density as the sum of the thermal density nth at 
the cryostat bath temperature and the excess density riox- 
Relevant quasiparticle distributions are strongly peaked 
at energy A, thus for the purpose of tracking energy dur- 
ing quasiparticle relaxation we treat all quasiparticles as 
having energy A. When we describe the superconductor 
as having a temperature Tg 7^ Tj,, we choose Tg to have 
the correct quasiparticle density. 
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FIG. 4: Model sweep over the value of the resistance area 
product for the overlayer trap barrier. Large resistance area 
product values represent the case of no overlayer trap, while 
low values represent a highly transparent trap barrier. The 
right vertical axis shows values of Tn, the left vertical shows 
the power loads from Eq|4]as equivalent /3 values. The solid 
red terms sum to equal 'Ideal Pn'- 'Ideal Pn' is -PN|Ts=Tb ^nd 
'Lost Pn' is the difference between 'Ideal Pn' and the actual 
Pn due to heating of the S layer. The terms with x markers 
are often represented as pPs. Decreasing the trap resistance 
area product causes the S layer temperature Ts to remain 
closer to the bath temperature Tb, reducing 'Lost Pn' and 
recombination, and resulting in lower base temperature Tn. 
Changes in Tn lead to changes in 'Ideal Pn', the electron- 
phonon coupling, and the power load due to the athermal 
behavior of trapped quasiparticles. The majority of the ben- 
efit of the overlayer trap is achieved with a resistance area 
product of 10 — 100 f2/im^, where Tn values near 75 mK are 
reached from Tb = 300 mK. 'Lost Pn' (and therefore /3) re- 
mains finite even for low resistance area product overlayer 
because the overlayer trap heats as well. The Andreev (hVh) 
power load is not shown; it is roughly constant and never 
greater than the recombination term. 



B. NIS Tunneling 



The form and magnitude of the two particle tunneling 
current depends on electron interference due to multiple 
scatterings, and therefore, on the geometry of the elec- 
trodes. We rely on theoretical forms of /2|n and /2|s 
from Hekking and Nazaro\n^ which are not specific to 
our junction geometries. The form of I2 in Eqs[9|and |10| 
is in rough agreement with our data, and we have an on- 
going investigation to improve our understanding of the 
precise form of I2 for our geometry.'^ We have excluded 
the contribution /2IS from the superconducting electrode 
in this work because 1) J21S -^2|n due to the thick- 
ness of the superconductor 2) the overlayer traps should 
further suppress multiple reflections and thus the mag- 
nitude of /21s and 3) the theory of /21s in Hekking and 
Nazarov'^ breaks down for biases QeVh ~ A, which are 
commonly used for cooling. In most cases Ii 3> /2, how- 
ever at low temperatures and low biases I2 > h- For 
the junctions described in this work I2 plays a small role; 
I2 is more important for junctions with lower resistance 
area products. 

The complete current-voltage (IV) relationship also in- 
cludes a resistive voltage due to the normal metal elec- 
trode 

V = l^b+/NS^^pad, (11) 

The power deposited in the normal metal by single 
particle tunneling is 

1 



Qe-^NS J- 

-h{E)]dE 



(qel4 - E)v{E)[fn{E - q,V^) (12) 



Refrigeration is possible because Pn is negative for biases 
such that (jeVb ~ A. Note that Pn is a function of both 
Tn and Tg (through /n and /g), and this dependence 
on Tg is the reason that superconductor heating directly 
impacts NIS performance. The two particle tunneling 
current deposits power l2Vh and Joule heating deposits 



power I^a^R-pad in the normal metal.E! 



The current in an NIS junction is made up of two parts The power deposited in the superconductor is 
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(9) 
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Pg = /iFb - Pn 



(13) 



where Ii is the single particle tunneling current, I2 is the 
two particle tunneling current, /2IN is the two particle 
contribution from the normal metal electrode in the sim- 
plest geometry of an infinite uniform junction, /2|g is 
the two particle contribution from the superconducting 
electrode, Qe is the electron charge, Pns is the tunneling 
resistance, TZyis is the product of the tunneling resistance 
and junction area, and Vb is the voltage difference across 
the junction. 



Both the quasiparticle thermal population, and the 
quasiparticles injected into the S layer by the NIS junc- 
tion are strongly peaked at A. We approximate the 
quasiparticle injection rate as PgA^^, which is justified 
because Tn ^ A/fcb and eVb < A in the regime of in- 
terest for NIS refrigerators. The power deposited in the 
superconductor per unit area Pg is calculated by sub- 
stituting the resistance area product 7?.ns of the tunnel 
junction in place of the resistance Pns- 



IV. QUASIPARTICLE TRAPPING 

A. NIS junction as a quasiparticle trap 

Power will flow across an NIS junction even when 
Vb = 0, if the normal metal temperature Tn and su- 
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perconductor temperature Tg are unequal. In NIS refrig- 
erator operation, the S layer is directly heated and there- 
fore is hotter than the O layer. As a result, power will 
flow from the S to the O layer, providing an additional 
mechanism for the superconductor to reach thermal equi- 
librium. The unbiased NIS junction between the S and 
O layers, combined with the O layer form a quasiparti- 
cle trap used to heatsink the superconductor. The power 
flow per unit area from the superconductor is given by 



Vtr 



ap 



9o^trap Jo 

A 

g2^trapiV(0) 



Eu{E)[fsiE) - h{E)]dE (14) 
[n{Ts) - n{To)] (15) 



where we have used the approximation that all quasipar- 
ticles have energy A to obtain Eq[T5| 

The tunneling lifetime of an electron in a thin metal 
film adjacent to a tunnel barrier is 



(16) 



where t is the film thickness and TZ^ is the resistance area 
product of the tunneling barrier. The tunneling lifetime 
grows with film thickness and also with the resistance 
area product of the tunnel junction. The lifetime for 
quasiparticles to trap from the S to the O, or to tunnel 
back from the O to the S layer is Tt — 104 ns, based 
upon thickness t « 500 nm, and a tunnel barrier with 
resistance area product 7?,t = 60 fl^m'^ . 



B. What happens to a trapped quasiparticle? 

A quasiparticle which tunnels from the S into the O 
layer becomes an excited electron with energy A. There 
are three processes available to this electron: 1) tunneling 
back into the superconductor with lifetime rtun — 104 ns, 
2) scattering with another electron with lifetime Tc_o|a = 
540 ns, and 3) scatter and create a phonon with lifetime 
''o-pIA = 25 ns. Once the excited electron has scattered 
with either a phonon or an electron, it will have energy 
less than A and be unable to tunnel back into the super- 
conductor, thus the name quasiparticle trap.l^ However, 
the most likely phonon to be created will have energy 
3 A / 4 which is well above the thermal distribution.li^ A 
3 A / 4 phonon has a finite probability of being absorbed in 
the N layer, and if this happens regularly it will severely 
degrade the benefit of quasiparticle traps. Appendices 
lA IHA 4l describe these mechanisms and the methods used 
to calculate the time constants. 

Wc estimate the probability of scattering to create a 
phonon Hp, to tunnel back to the superconductor Iltunj 
or to undergo electron-electron scattering Ho as 



n 



(p)/(tun)/(c) 



'(c-p|A)/(tun)/(o-o|A) 



-p|A 



(17) 



-e|A 



For A = 190 fJ,eV, Hp = 71% of the quasiparticles re- 
moved from the S layer will generate athermal phonons, 
Iltun = 24% tunnel back to the superconductor, and only 
He — 5% are thermalized by electron-electron scattering. 
The most likely energy for these athermal phonons is 
E = 3A/4, so - I Hp = 53% of the energy which tunnels 
from to O becomes athermal phonons. Appendix |C 4| 
describes a ray tracing model which predicts the frac- 
tion of these athermal phonons that deposit their energy 
in the other layers. The majority of athermal phonons 
will be reabsorbed before leaving the O layer, making 
the overlayer traps more effective than these probabili- 
ties suggest on their own. 



V. EXPERIMENTAL METHODS AND RESULTS 

SINIS refrigerator devices were fabricated on a Si wafer 
with ~150 nm of thermally grown Si02. The N layer 
is ~20 nm of sputter deposited Al with ^ 4000 ppmaPll 
Mn. The N layer and all subsequent layers are patterned 
with standard photolithography techniques. The metal 
layers are etched with an acid etch. A layer of Si02 is 
deposited by plasma enhanced chemical vapor deposition 
and etched with a plasma etch. The vias etched in the 
Si02 define the location of the tunnel junctions. The 
wafer is returned to the deposition system where it is ion 
milled to remove the native oxide from the AlMn, then 
exposed to a linearly increasing pressure of 99.999% pure 
O2 gas for 510 seconds, reaching a maximum pressure 
of 6.0 Torr (800 Pa) to form the tunnel barrier oxide. 
Then, ^500 nm of Al is deposited by sputter deposition 
to form the S layer. The overlayer trap tunnel barrier 
is formed by exposure to 32 mTorr O2 (4.3 Pa) for 30 s, 
and ^500 nm of AlMn is deposited to form the O layer. 

One of these devices (inset in Fig [5]) was cooled in a 
copper box attached to the cold stage of an adiabatic de- 
magnetization refrigerator. We measured current-voltage 
{IV) curves of the SINIS device (and the independent 
thermometer junctions) at many bath temperatures Tb. 
The bias of the refrigerators (independent thermometers) 
was set with a computer-controlled, battery-powered 
voltage source connected through a 100 kfi (10 Mft) bias 
resistor. The voltage across the device was amplified with 
an operational amplifier and the output of the amplifier 
was measured by a digital multimeter. The current was 
calculated as the source voltage minus the device voltage 
divided by the bias resistance. 

The primary experimental results we report come from 
thermometry based on the refrigerator junctions. Each 
data point (/ and V) is an independent measurement 
of the normal metal temperature which depends on 
both bath temperature Tb and refrigerator bias V . The 
temperature Tn at each point is uniquely determined by 
comparison to isothermal theory curves (Eqs[7]and[ll]), 
as shown in Fig [5j Temperatures determined from NIS 
junctions in this way are effective temperatures because 
Fermi occupation functions are assumed in the junction 
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FIG. 5: Theoretical isothermal NTS current-voltage curves 
from Eqs [7| and 1 1 1 1 shown as solid lines. Measured NTS refrig- 
erator IV curves taken at 100 mK and 300 mK bath temper- 
atures shown as circles. For each data point, the temperature 
Tn is uniquely determined by the theory curve upon which 
it falls. The optimal bias in the data for each bath temper- 
ature is shown as an x marker. The inset shows a picture 
of the device. The two large junctions (vertical rectangles) 
are the refrigerator junctions. The two smaller junctions are 
independent thermometer junctions. 
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FIG. 6: Temperature Tn determined from the refrigera- 
tor junctions vs NIS refrigerator voltage bias as solid lines 
for many bath temperatures Tb. The vertical dashed line at 
20 /iV is used to determine uncooled NIS data points, for use 
in calibrating the conversion between Vf and Tn. The optimal 
bias points are shown as x, and indicate lowest temperature 
Tn achieved at each bath temperature Tb. 

electrodes. For a thermal electron distribution this effec- 
tive temperature will match the actual temperature. For 
athermal electron distributions, the NIS effective temper- 
ature does not provide a complete description. Figure [6] 
shows effective temperature Tn vs refrigerator bias for 
many bath temperatures. 

For each bath temperature, we determine the bias 
point at which the minimum value of Tn is reached. We 



call this the optimal bias. The value of Tn at the optimal 
bias is shown vs bath temperature in Fig [7] Figure [7] also 
shows the value of Tn at a bias voltage of 20 fj,V. For 
bath temperatures above ~ 250 mK the temperature Tn 
determined at a bias voltage of 20 fiV should equal the 
bath temperature because the NIS junctions have little 
thermal effect at low biases. We determine A by a least 
squares minimization between the temperature deduced 
from IV curves at low bias (20 fiV) and the temperature 
measured by the cryostat thermometer. The N layer re- 
sistance i?pad is calculated as in Table |Tj the tunneling 
resistance Rt is the asymptotic resistance at high bias 
and high temperature minus i?pad, and 7 is chosen to 
match the maximum differential resistance at 70 mK. 

The data from the refrigerator IVs indicate signifi- 
cantly improved cooling compared to previous large area 
(> 100 /im^) NIS refrigerators.'^ Assuming a thermal 
electron distribution in the normal metal, we deduce cool- 
ing from 300 mK to 82 mK and from 100 mK to 26 mK. 

VI. INDEPENDENT THERMOMETERS AND 
ATHERMAL ELECTRON DISTRIBUTIONS 

Current-voltage curves from the independent ther- 
mometers can be used to obtain temperature values us- 
ing a procedure similar to that described in the previous 
section. As shown in Fig 7, the thermometer data sug- 
gests hotter electron temperatures than those deduced 
from the refrigerators. The thermometer junctions were 
fabricated using a double oxidation technique and their 
resistance area product was 145,0 00 O /im^, over 100 
times greater than the refrigerators This high resis- 
tance makes the thermometers thermally neutral, mean- 
ing they neither heat nor cool the normal metal. We next 
consider thermal mechanisms for a temperature gradient 
between the thermometers and refrigerators. 

We expect a small temperature difference between the 
refrigerator and thermometer junctions due to the fi- 
nite thermal conductivity of the normal electrode and 
the presence of power loads within it. We calculate 
the expected difference by solving the heat equation in 
the N layer including electron-phonon coupling and stray 
power: 

^ - S(T^ - T6) - Texccss/t/N (18) 

We let X vary from to 10 microns with a; = corre- 
sponding to the edge of the refrigerator junction, a; = 10 
corresponding to the end of the AlMn, and x = 4 to 
9 corresponding to the thermometer junction. We fix 
dT/dx at x = at a; = 10 and vary the refrigerator tem- 
perature at a; = until the thermometer temperature 
(evaluated at x — 8) matches the measured tempera- 
ture. The variable kn is the thermal conductivity of the 
normal layer and is obtained from Eq |A2[ The stray 
power Pexcoss = 0-05 pW is obtained from thermometer 
junction temperature measurements with the refrigerator 
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junctions at zero bias. Finally, L^n is the N layer volume. 
The observed thermometer temperatures and the refrig- 
erator temperatures calculated from these thermometer 
temperatures are shown in Fig [7| The gradient varies 
from 0.02 mK to 7 mK for bath temperatures between 
100 mK and 500 mK. This gradient accounts for over 
half of the observed difference at a bath temperature of 
500 mK but does not explain the observed difference near 
100 mK. We hypothesize that the remaining differences 
in temperature result from an athermal electron distribu- 
tion created by the refrigerators that recovers to a ther- 
mal distribution over the distance to the thermometers. 

NIS refrigerators generate an athermal electron distri- 
bution when the tunneling rate is faster than inelastic 
scattering in the normal electrode. The tunneling time 
for electrons above the gap edge is insensitive to temper- 
ature, and the rates for inelastic electron-electron and 
electron-phonon scattering both decrease with tempera- 
ture. As a result, athermal effects in NIS junctions ap- 
pear near or below ^ 100 mK.I^ The athermal electron 
distribution created by the refrigerators lacks the high 
energy excitations ordinarily found in the tips of a Fermi 
distribution. However, its total excitation energy will ex- 
ceed that of the thermal distribution that results in the 
same tunneling current. 

Over the distance between the refrigerator and ther- 
mometer junctions, the athermal distribution recovers to 
a thermal one through inelastic scattering and the tips 
of the Fermi distribution are repopulated. Consequently, 
additional current flows through the thermometers and 
a higher device temperature is deduced even though the 
thermal electron distribution at the thermometers and 
the athermal distribution at the refrigerators contain the 
same excitation energy (neglecting the small gradients 
calculated above). We observe that the temperature de- 
duced from the thermometers is largely independent of 
their bias point, supporting the notion of a thermal elec- 
tron distribution at the thermometer junctions. Further, 
athermal electron distributions in metals with similar dif- 
fusion constants to our normal layer have been observed 
to thermalize over distances as short as 2.5 

Our results illustrate the complexities of junction ther- 
mometry. Experiments that use the same junctions as 
both refrigerators and thermometers are susceptible to 
athermal effects. Similarly, if independent thermome- 
ter junctions are close enough to the refrigerators that 
they sample the same athermal distribution, the inferred 
temperature reduction may be exaggerated. We have 
provided temperature values deduced from both refriger- 
ator and distant thermometer junctions. While the ther- 
mometers likely provide a truer measurement of temper- 
ature, the results from the refrigerators still have value 
as a measure of the distortion of the electron distribution 
in our devices. 

We note that our equilibrium thermal model predicts 
temperatures close to the values deduced from the re- 
frigerator IVs, rather than the independent thermome- 
ters. This outcome is not surprising. While the model 



finds the temperature where the power loads balance, it 
can also be thought of as finding the junction current at 
thermal balance. In the presence of an athermal elec- 
tronic distribution, current remains a good predictor of 
key power loads such as the Joule term and the heating 
of the superconductor. Hence, the model correctly pre- 
dicts refrigerator currents together with the temperature 
that produces these currents assuming a thermal distri- 
bution. Since the refrigerator temperatures are deduced 
from measured currents assuming a thermal distribution, 
model and refrigerator measurements give similar tem- 
perature values. 

VII. ANALYSIS AND COMPARISON TO 
MODEL 

Since knowledge of the precise details of the electron 
distribution under refrigeration is lacking, for the pur- 
poses of this section we henceforth assume a Fermi dis- 
tribution in the normal metal and continue our analysis of 
the results from that perspective. We calculate /3 vs bath 
temperature and compare to predictions of the thermal 
model, results are shown in Fig|8] In order to determine 
the value of /3 in the device, we calculate the excess power 
load required to explain the measured value of Tn at the 
optimal bias at each bath temperature; 

Pcxccss = - 2PN(i?T/2, A, Vo/2, Tn, Tb) ^^^^ 
— loRpad — Ujs!T,{T^ — T§)) — I2V0 

where Rt is the tunneling resistance of the two refriger- 
ator junctions in series, the factors of 2 account for the 
existence of two junctions, Pn is given in Eq[T2] Vo and 
lo are the optimal bias voltage and current for a given 
bath temperature, Tn is the N layer electron tempera- 
ture at the optimal bias, i?pad is the resistance of the N 
layer current path, Ujq is the N layer volume, and S is 
the electron phonon coupling constant. The value of /3 is 
calculated by 



/oT4-2PN(i?T/2,A,T/o/2,TN,Tb)' ' ' 

where IqVq is the dissipated power, and the denominator 
is the total power deposited in the S layer. The mea- 
sured value of /3 is therefore a complicated combination 
of measurement and analysis that depends on the param- 
eters A, S, 7, [/n, Rt and i?pad, and the accuracy of the 
underlying theory. We generate error bars on /? using 
estimates of the uncertainty of each parameter propa- 
gated through the calculation. The uncertainties used 
to generate the error bars are A : ±0.25%, S : ±20%, 7 : 
±10%,J7n : ±5%,andPpad : ±10%. Thesumof Ppad and 
-Rt was held constant. Variation in A affects the deter- 
mination of Tn and directly enters the calculations of Pn. 
At higher temperatures uncertainty in S determines the 
bulk of the error bars; at lower temperatures uncertainty 
in A dominates. The range of 300 mK and below is the 
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FIG. 7: Temperature of the normal metal Tn determined 
from refrigerator junctions at optimal bias voltage as o vs 
bath temperature Tb. Temperature Tn from the indepen- 
dent thermometer junctions as T, and the temperature at 
the refrigerator junction calculated from thermometer junc- 
tion temperature as A (see Sec VI I. Uncooled temperature 
Tn at Vf = 20 fj,Y as 0- The solid blue line is the predicted 
temperature Tn at optimal bias and the solid grey line is 
Tn = Tb, which represents zero heating or cooling. The in- 
set shows the difference between the model and measured Tn 
at optimal bias as o, and a subset of the difference between 
the temperature Tn at Vp = 20 and Tb as <). Uncooled 
temperature points at bath temperature below 250 mK are 
excluded because they are extremely sensitive to the form of 
the two particle tunneling current, which is not fully under- 
stood. The theory and data for Tn at optimal bias are in 
very good agreement. These results are an improvement over 
previous cooling with large area NIS junction refrigerators!^- 



most interesting range for Al-based NIS refrigerators. In 
this range, the measured value of /3 is roughly 1.5 times 
the predicted value, which is better agreement than past 
NIS modeling has achieved. For higher temperatures, 
the agreement is not as good, suggesting that improve- 
ment could be made by focusing on features of the model 
whose relative importance increases with temperature, 
such as quasiparticle recombination and electron-phonon 
coupling. Underwood et aP^l observe that the exponent 
in the electron phonon coupling depends on temperature 
for temperatures above ^ 300 mK; including this effect 
may improve agreement at higher temperatures. 



VIII. FUTURE REFRIGERATOR DESIGNS 
BASED ON THERMAL MODEL 

By exploring the design parameter space with the ther- 
mal model, we have developed an NIS design that is pre- 
dicted to cool from 100 mK to 6.5 mK. The overlayer 
trap resistance TZq is decreased to 1 rtfim?, the overlayer 
thickness to is increased to 10 /im, the Dynes parameter 
7 is reduced to the lowest reported value^-of 38000"^, and 
the junction length parallel to current flow Xj is reduced 



FIG. 8: Measured values of /3 vs Tb as circles and model pre- 
dictions as solid line. The inset shows the ratio between the 
measured and predicted value vs Tb. Previous predictions of 
/3 or equivalent parameters have varied from measured values 
by factors on order ~ 3 — 10, so a factor of 1.5 is considered 
good. The model predicts the cooling to within a few mK at 
low temperatures, which is useful for designing improved NIS 
refrigerators. 



to 3 fim. 

The same changes that we suggest to improve cooling 
from 100 mK will also improve cooling from higher tem- 
peratures with lower resistance area product refrigerator 
junctions. A device with the changes above, increased 
N layer thickness of — 30 nm, and decreased the re- 
sistance area product of 7?.ns=150 Q^m^ is predicted to 
cool from 300 mK to 48 mK. If the overlayer thickness 
is left unchanged at 500 nm, the device should cool to 
64 mK. Junctions with this low resistance area product 
are difficult to fabricate with high yield, so it is not a 
simple step to achieve these performance levels. These 
predictions demonstrate the capability of the model to 
explore a large design space. Refrigerators based on these 
designs will generate electron distributions that vary fur- 
ther from the Fermi distribution than the refrigerator 
discussed here. 



IX. CONCLUSIONS 

We have described a comprehensive NIS refrigerator 
thermal model that includes differential equations for 
three temperatures vs position, and accounts for the 
athermal behavior of high energy excitations. The pre- 
dictions of this zero free parameter model agree well with 
experimental results from an NIS refrigerator with over- 
layer traps. The agreement for the quasiparticle return 
parameter P are within a factor of 1.5 over the temper- 
ature range of primary interest. This refrigerator cools 
more than previous large area (> 100/im^) NIS refrigera- 
tors. From a bath temperature of 300 mK, the refrigera- 
tor junction current is consistent with cooling to 82 mK. 
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However, the refrigerators likely generate an athermal 
electron distribution, and as a result the cooling can- 
not be described simply. Thermometer junctions which 
probe the temperature ~ 8 fim away from the athermal 
distribution measure a temperature 114 mK. 

This work expands upon previous modeling efforts by 
including more coupled temperatures systems (the over- 
layer electrons and the metal layer phonons) and more 
athermal effects. The model could be further improved 
by including more detail in the form of the electron- 
phonon coupling, and by calculating the electron distri- 
bution in the normal metal electrode. We plan to perform 
phonon cooling experiments that will provide a direct 
measurement of the usable temperature reduction from 
similar refrigerator junctions. 
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where L = 2.45 x 10 * W ft/K'^ is the Lorenz number.'i^ 



3. Electron-phonon interaction in normal metals 

Electron and phonon populations occupying the same 
volume of a normal metal interact via the electric field 
from phonons displacing lattice ions .^^^E^ This interaction 
leads to a power per unit volume flowing from the phonon 
system to the electron system of the form 

7'p-e = S(T?-T;"), (A3) 

where Tp is the phonon temperature. To is the elec- 
tron temperature, E is the electron-phonon coupling con- 
stant, and n is an exponent whose value can vary from 
4 to 6 depending on the properties of the metal film. 
Previous work on NTS junctions with AlMn base elec- 
trodes has used n = 5!^ However, recent measurements 
by Underwood^ and Taskinen^'^ show that n = 6 is more 
accurate for AlMn films in the temperature range of in- 
terest. 



a. Electron-phonon interaction: Electron scattering time 



Appendix A: Electron excitation behavior in a 
normal metal 



To calculate the temperature of the S layer, the O 
layer, and the effect of athermal excitations in both lay- 
ers, we must understand the heat transport and loss 
mechanisms in both a superconductor and a normal 
metal. In particular, we want to understand what hap- 
pens when a quasiparticle is trapped in the O layer. 



1. Electron diffusion in the normal metal 

The numerous elastic collisions off impurities, grain 
boundaries, and surfaces cause electrons to diffuse, rather 
than move ballistically, through a normal metal. The dif- 
fusion constant for electrons in a normal metal i^^^ 



1 



(Al) 



where /x is the elastic mean free path and vp = 2.02 x 
10^ m/s is the Fermi velocity for Al. From this we calcu- 
late the elastic mean free paths in the layers of the NIS 
refrigerator to be = 6 nm, Iq = 9 nm and /g — 340 nm. 



2. Thermal transport in normal metal electrons 

The Wiedemann-Franz law relates the thermal conduc- 
tivity of the electron system to the resistivity by 



Px 



(A2) 



UUorrfl^ extended theoretical calculations for Pp_o to 
find To-p, the characteristic time for an electron to down 
scatter with the creation a phonon. In the case where 
n—5 in Eq|A3| 



'^c— p,n— 5 

where /o=r(5)C(5) w 25. 



3/ofc^A^(0) 



(A4) 



To obtain rc_p for the n — 6 case we use an energy scal- 



ing argument to modify Eq |A4[ A method to estimate the 
electron lifetime for phonon scattering is r^-p « Ce/Qc-p 
where Co = is the electron heat capacitjff^ per unit 
volume and 



^c-p = dV,-p/dT, = nET^- 



(A5) 



is the electron-phonon thermal conductivity. Thus 
To_p ~ 7n(j/(riSr"~^). This method underestimates 
To_p because it represents an average of To_p over ther- 
mal occupation, and the energy dependence of r^-p is 
strong. However, this method accurately predicts the 
energy scaling of r^-p. Therefore we modify Eq X4 for 
n = 6 by replacing Y.E^ with SS^/fcb. The result is 



3IoklN{0) 



-p,n= 



(A6) 



We introduce the notation r, 

T, 



-p|A = 25 ns, to represent 



„_p evaluated at energy A. 
An electron scattering with emission of a phonorfi^l 
with initial energy A will have average final energy A/4 
with the emission of a phonon of energy 3A/4.'321 Since 
this phonon's energy is well above the thermal distribu- 
tion at 300 mK its athermal behavior must be accounted 
for, as described in Appendix |C 4[ 
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Electron-phonon interaction: Phonon scattering time ing times are given by 



The energy dependence of the hfetime Tp_o for a 
phonon of energy E to lose energy by scattering and pro- 
moting an electron is weaker than for Te_p, so the ap- 
proximation Tp_e ss Cp/^p_o is more appropriate. Using 
the Debye result for phonon heat capacity, 



234iVAndA:i^ 



(A7) 



^c—e, clean, 3D 
To-e, dirty, 2D 



qlpE\n{KT^/E) 



(A8) 
(A9) 



where Ti « 10^^ K. When an electron with initial energy 
A scatters, the average final energy is A/3 with the cre- 
ation of two new excitations, an electron and a hole, each 
with energy A/S.SSiSl^ 



4. Electron-electron scattering 

The electron-electron scattering rate t~^^ is strongly 
dependent on the effective dimensionality of the sample, 
and the degree of disorder present. Therefore, we will 
first look at the criteria for the 2D/3D limits and the 
clean/dirty limits. 



a. 2D vs 3D scattering 



Appendix B: Excitation behavior in a 
superconductor 

Having looked at the mechanisms of thcrmalization 
that take place in a normal metal, we now examine their 
analogs in a superconductor. Quasiparticles undergo dif- 
fusion by many elastic collisions, much like electrons, but 
inelastic scattering is rare for quasiparticles with energy 
~ A. Instead, the most common method for quasiparti- 
cle relaxation is recombination, whereby two quasiparti- 
cles combine to form one Cooper pair and one athermal 
phonon with energy 2A. 



The material is in the 2D limit if the excitation energy 
of the electron is less than the uncertainty in the energy 
associated with the transit time across the thinnest di- 
mension, that is Ttiansit < h/E. Otherwise it is in the 
3D limit. All of the films here are thicker than the elas- 
tic mean free path so the transit time is Ttransit = t'^ /D 
where D is the diffusion constant. Therefore Ttiansit-N = 
1.1 X 10-13 s and Ttransit-O = 4.0 X 10"" s. The ex- 
citation energy of most interest for scattering is A, and 
h/A = 3.5 X 10-12 s.ESIWe determine that the O layer 
is in the 3D limit and the N layer is in the 2D limit for 
electrons with energy A, such as trapped quasiparticles. 



b. Electron- electron scattering: clean vs dirty scattering 

The relevant condition for clean vs dirty electron- 
electron scattering limits is different in the 2D and 3D 
cases. In 3D, dirty scattering dominates when E < 
E-p/ik-plY" and in 2D dirty scattering dominates when 
E < {hvF/l){n/{kFt)). The N layer is in the dirty limit, 
and the O layer is in the clean limit. 



c. Electron- electron scattering time: r^-e 

The O layer films look 3D and clean to electrons of 
energy A, and the TV layer films look 2D and dirty to the 
same electrons. Therefore the electron-electron scatter- 



1. Quasiparticle diffusion 

We treat all quasiparticles as having equal energy A 
with one exception. This exception is the quasiparticle 
group velocity, which is strongly energy dependent and 
approaches zero as E approaches A. The group velocity 
ws {E) for a quasiparticle excitation of energy E i^MMMI 



vs{E) 
(«s(Ts)> 



VF^l~{A/E) 
2iV(0) f 
<m Jo 



(Bl) 

vs{E)UE)ya{E). (B2) 



The diffusion constant for quasiparticles in a supercon- 
ductor is modified compared to the normal state by this 
group velocity, such that Z3s = -Ds-normai(i's)/'f^F- How- 
ever, the quasiparticle distribution most relevant to NIS 
refrigeration is not thermal. 

An NIS refrigerator junction injects an athermal dis- 
tribution of quasiparticles into the S layer. The diffusion 
constant for these injected quasiparticles is 



-^S— norm; 

rs(i4) 

h{E - q.Vx,) - 2h{E)]dE 



VY (B3) 



where Fs is the quasiparticle injection rate, which can 
be calculated with this same integral with ws/vp — >■ 1. 
For an NIS refrigerator cooling from 300 mK to 100 mK, 
£'s-normai/£'si = 20 and D^-i = IQi" /iTO^/s are typical 
values. 
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2. Quasiparticle-phonon scattering 

Like an excited electron in a normal metal, a quasi- 
particle can relax to a lower energy state by emitting 
a phonon. The lifetime for this is Tqp_p. High energy 
{E >> A) quasiparticles relax at a rate almost identical 
to high energy electrons in a normal metal. For energies 
close to A, quasiparticle-phonon scattering is rare. In Al 
at a temperature of 250 mK and an excitation energy of 
1.2A, the scattering time Tqp_p is many microseconds. 
Therefore, Tqp_p will be much longer than other relevant 
time constants and quasiparticle-phonon scattering is not 
considered for the NIS refrigerators discussed her^i^. 
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3. Quasiparticle-quasiparticle recombination 

Two quasiparticles can recombine to form a Cooper 
pair and emit a phonon with energy equal to the sum 
of their excitation energies. This will result in athermal 
phonons with energy 2 A, whose behavior is described in 
Appendix |C 4| Since recombination is a two-body pro- 
cess, the recombination rate scales as the number of pair- 
ing possibilities 



(B4) 

The coefficient Fr 
(B5) 



dt 

where n is the quasiparticle density, 
depends on material properties, 

_ / ^ ^,3 4 
" ^k^J A^s(0)ATo_qp 

where To-qp is a material specific parameter, whose value 
we take to be 100 ns in the S layer. In order to compare 
the recombination rate to other time constants, it is use- 
ful to define Tqp-qp = 1/ (Frtt.). For thermal quasiparticle 
populations at 300 mK and 400 mK, Tqp_qp = 7.7 ^s and 
1.0 fis respectively. Recent measurements of quasiparti- 
cle lifetime in Al by Barends et aP^ provide a check to 
these calculations. At 210 mK, they measure a lifetime 
of about 180 /iS and we calculate rqp_qp = 210 /xs, which 
is in very good agreement. 

In the thermal model, we separate the total quasipar- 
ticle density into two parts, the thermal density nth eval- 
uated at temperature Tb and the excess density riox- The 
excess density is an effective increase in the superconduc- 
tor quasiparticle system temperature due to the power 
deposited by an NIS junction. The rate of excess quasi- 
particle recombination is given by 



dt 



(B6) 



Phonon-quasiparticle scattering or 
breaking" 



pair 



The excitation of electrons by phonons has already 
been discussed in a normal metal. A similar process 



FIG. 9: The propagation of a phonon starting at position Zi 
in the O layer and propagating at angle ip. The probability 
amplitude of the phonon is reduced as it propagates through 
each layer depending on the time constant r for the phonon 
to interact in that layer. Additionally there is an escape prob- 
ability (ri) at the interface between the bath and layer A''. By 
following the phonon until the amplitude is reduced to near 
zero, we calculate the probability that a given phonon either 
escapes to the bath or is absorbed in the O, S or N layer. 



can occur in a superconductor where a phonon of energy 
2A or greater destroys a Cooper pair and produces two 
quasiparticles. When two quasiparticles recombine, the 
emitted phonon has sufficient energy to break another 
Cooper pair. The pair breaking time for a 2 A phonon 
is given by KaplarP^ and is equal to a material specific 
parameter Tp_pb = 230 ps for Al. This number is only 
weakly dependent on temperature and phonon energy. 
For example, the rate is only 20% faster for a 3A phonon 
compared to a 2A phonon. 



Appendix C: Phonon escape 

Ideally, phonons generated by quasiparticle recombina- 
tion or electron-phonon scattering escape into the Si sub- 
strate, but may instead interact with electrons or quasi- 
particles before leaving the Al and AlMn layers that make 
up the NIS refrigerator. 

While the electrons in the N, S, and O layers are iso- 
lated from each other by tunnel barriers, the phonons are 
not similarly isolated. All three layers are Al with only 
small quantities 4000 ppma) of Mn in the N and O 
layers, so phonons should freely pass between these lay- 
ers. However, the metal layers making up the NIS refrig- 
erator are deposited on ~ 150 nm of Si02 on a Si wafer 
300 /im thick, and because of this material mismatch the 
phonons will not always escape into the Si wafer. The Si 
wafer is treated as an ideal thermal bath at temperature 
Tb. We model the behavior of thermal phonons using an 
acoustic mismatch boundary resistance, and the ather- 
mal phonons using a ray tracing model. 
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1. Thermal phonons: acoustic mismatch 

Phonon thermal resistance at the interface between 
two different materials arises from reflection due to dif- 
ferences in the speed of sound and is known as acoustic 
mismatch. The power flow per unit area from phonons in 
the combined metal layers at temperature Tp to a Si02 
layer at temperature Tb is 

Pamm = - T^) (CI) 

where ^ = 360 pW/ (/im^ K^) for Al on SiOsP 

2. Thermal phonons: thermal conductivity 

The thermal conductivity Kp of the phonons in the 
combined metals layers can be calculated as CpZp(s)/3 
where Cp is the Debye model heat capacity for phonons, 
Ip is the mean free path, and (s) is the average phonon 
speed. Boundary scattering limits the phonon mean free 
path at low temperatures, so we approximate Ip = tg + 
to +tN- Therefore 

(iN +to+ <s)(s>234iVAndA:b , T 3 

where the average speed of sound (s) = {s\ + 2st)/3 = 
4.4 X 10'^ m/s is a weighted average of the longitudi- 
nal si = 6.7 X 10'^ m/s and transverse St — 3.3 x 
10"^ m/s values for Al. One could use a different value 
of (s) for phonons originating in the superconductor 
(primarily longitudinal) or the normal metal (primarily 
transverse) ,23^1 the choice of weighting for (s) has 
only a small impact on the results of the thermal model. 
The extreme cases of (s) ~ st and (s) = si result in val- 
ues of Tn that differ by ~ 0.5 mK, for the parameters in 
Table HI 



3. Phonon-phonon scattering 

Phonons can relax by splitting into two phonons of 
lower energy, but this anharmonic decay process is much 
slower than other relaxation pathways. Mari^^H states 
that the lifetime for anharmonic phonon decay is of order 
Tp_p fa 10 Y>s{khQu/E)^ , which is ~10 ms for phonons 
of energy 2A in Al. There is experimental evidence for 
phonon mean free paths of order 1 mm in superconduct- 
ing materials, consistent with long lifetime for anhar- 
monic decay.l^ 

4. Athermal phonons 

Both quasiparticle recombination and electron-phonon 
scattering by trapped quasiparticles will create athermal 
phonons. These phonons will interact with electrons or 



escape before relaxing via anharmonic decay, and there- 
fore their behavior will not be well described by equilib- 
rium physics. In this section, we describe a ray tracing 
model shown in Fig[9]for calculating the behavior of these 
phonons. ■^■^ With this model, we calculate the fraction of 
these phonons which are absorbed in each metal layer, 
and the fraction which escape to the bath. The frac- 
tion of energy deposited in the N layer is of particular 
interest because of its direct impact on NIS refrigerator 
performance. 

Consider a phonon with energy 3 A /4 created at the top 
of the O layer traveling directly downward with speed of 
sound (s) . The time to propagate from the top of the O 
layer to the OS interface is to / (s) • The probability that 
the phonon scatters with an electron in the O layer in this 
time is 1 — e~*°/(^*^'^p— =13^/''). If the phonon enters the S 
layer, the probability that it leaves without interacting 
is 1 since its energy is below 2A. The probability that 
it interacts in the N layer before reaching the iV-Bath 
interface is 1 - e"*"/(<''>^p-=i3A/4). The probability that 
a phonon incident on the A^-Bath interface escapes to 
the bath is (77). Therefore, the probability that a phonon 
created at the top of the O layer escapes to the bath on its 
first attempt is e-*o/«''>-^p-=i3A/4)e-*N/((s)Tp_<,|3A/4)j^, xhe 

phonon escape probability (77) = (771 + 2rit)/3 = 0.71 is 
averaged over longitudinal and transverse phonons, and 
anglcli^l The value of the escape probability has a very 
weak effect on the results of the thermal model because 
most phonons are re-absorbed before they have even a 
single escape attempt. 

A phonon is created in either the O layer as shown, 
or in the S layer at position Zj and with direction given 
by angle (p relative to the junction. We assume specular 
reflections. The phonon is given an initial amplitude of 
1, which is reduced by a factor e"*/*-"^^"^ cosi^) jj-^ order to 
propagate the phonon to the next interface, where t is the 
distance to the interface and t is the time constant for 
the phonon to interact in the current layer. When the 
phonon reaches the interface between the N layer and 
the bath, the amplitude is reduced by a factor (1 — (77)) 
where (77) is the transmission probability for a phonon 
incident on an Al-Si02 interface. The probability of a 
phonon being absorbed in each layer or escaping to the 
bath is then calculated as the sum of all the amplitude 
reductions which occur in that layer, averaging over all 
possible initial values for z; and ip. 

The typical time for an electron of energy A to scatter 
with a phonon is T^-p\A ^ 25 ns. The average distance 
an electron will diffuse in this time is ^y Dqt^-pIa = 
13.3 fj,m, which is much greater than the overlayer trap 
thickness of 0.5 iim. Therefore all values of Zi are equally 
likely, although it would be necessary to restrict zj for 
thicker O layers or greater A because Tq-p\a falls as A 
increases. All values of (p are also equally likely. 

The variables that result from this calculation are 
called ^E-L where E is replaced with the energy of the 
phonon and L is replaced with the layer it is associated 
with. For example A2A-N is the probability that a 2A 
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phonon from recombination is absorbed in the N layer, from quasiparticle trapping is reabsorbed in the O layer, 
while ^3A/4-o is the probability that a 3/4A phonon 
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Input Parameters 


Parameter 


Input Value 




Description 


TZns 


1200 D,^im'^ 




junction resistance area product 


j4ns 


224 fim'^ 




area of a single refrigerator junction 




7 fj,m 




length of junction parallel to current flow 


^cdgc 


1 fj,m 




fabrication tolerance, contributes excess A*' volume 


pN 


0.0663 n^m 




resistivity of layer measured at 4 K 


tN 


17.5 nm 




thickness pf A*' layer 


^^sidc — trap 


3 /xm 




distance from junction to side-trap 


-^N — extra 


67.5 ^m^ 




parameter to account for excess A'' layer area 


E 


2.3 nW/i/im'^ 
60 Qfirri^ 




electron-phonon coupling strength 


7^o 




resistance area product of the overlayer traps 


Po 


0.0443 n/i"^ 




resistivity of the O layer measured at 4 K 


to 


500 nm 




thickness of the O layer 


Xcnd 


1000 /im 




length of overlayer trap 


A 


190 fieV 




energy gap in Al 


T-O-qp 


100 ns 




time constant related to quasiparticle recombination rate 




6000 




Dynes parameter, adds finite sub-gap states 


PS — normal 


0.00117 n^m 




resistivity of the S layer in the normal state, measured at 4 K 


ts 


500 nm 




thickness of the 5* layer 


Tb 


300 mK 




bath temperature provided by the cryostat 




360 pW/ifim'^ 




acoustic mismatch coefficient 




0.71 




athermal phonon escape probability 


(s) 


4.4x10^ /im/s 




average phonon speed in Al 


N{0) 


1.45 X 10^* 1/(J fim-^) 


two-spin density of states in Al 


Ep 


11.63 eV 




Fermi energy in Al 


Qd 


428 K 




Debye temperature in Al 


p 

J excess 


pW 




excess power added to the model 



Calculated 
Parameter 



Parameters 



Calculated Value 



Method of Calculation 



-Rns 

-^side — trap 
.^N— fridge 

C/n 

-fi'pad 



D. 



S — normal 



Ds^ 

«N 

Kp 



''"p-c|3A/4 

''"c — p| A 

^c—c\A 

Ttun-OS 

T"tun-SN 

^3A/4-0 

^3A/4-N 

A2A-O 

A2A-S 

Hp 

Iltun 
He 

g 

P'sidc — trap 



10' 
10" 



5.38 n 
10 fim 

323 fim^ 
6.83 pm^ 
0.59 n 
32 

2.3 X 10" fim^ 
9 X 10'-* fim^/s 
50 fj.m^/s 
1.66 X 10' 
3.70 X 
9.83 X 
9.2 ps 
66 ps 
38 ns 
540 ns 
111 ns 
2.2 /is 
0.86 
0.02 
0.33 
0.38 
0.17 
0.71 
0.24 
0.05 



7?.nsMns 

■^j ~t- f^sidctrap 

{xj + 2.5xcdgG)(yj + 2a;edgc), area of the A'' layer 

(^N + ^N-extra)tN 

pN(a;j/2-|- 1.5a;edge)/(tNj/j) 



^ W/(K Atm) 
W/(K nm) 
" W/(K fim) 



Eq 
Eq 
Eq 
Eq 
Eq 
Eq 
Eq 
Eq 
Eq 
Eq 
Eq 
Eq 



Al 



m 



-K7 



with Ps- normal 

with Ds- 



- normal 7 



Tn and Ts- 



see Appendix 
see Appendix 
see Appendix 
see Appendix 
see Ap pendix 
EqflTl 



with po and T = 300 mK 
with pN and T = 100 mK 
with T = 300 mK 
with E = 2A 
with £ = 3A/4 
with E = A 
with E = A 
with ts and TZo 
with ts and 7?.tn 



C4 



EqT7 

Eqp 

if X < Xj then g{x) — 1, else g{x) — 

if a: > Xsidc-trap then g{x) = 1, else g{x) = 



TABLE I: Numerical values and method used to calculate the parameters in the comprehensive thermal model. The input 
parameters are based upon measurements of the properties of the refrigerator described in this paper. The A'^ layer area depends 
on the aspect ratio of the tunnel junction due to fabrication tolerances which require a finite amount of AlMn to border the 
via defined junction. 



